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Abstract 

The production and propagation of light antimatter nuclei p, d, t, 3 He, ^He has been 
calculated using inclusive p production cross sections from a new data analysis, and 
coalescence models for the production of composite particles. Particles were propagated 
using recently proven phenomenological approaches. The atmospheric secondary flux is 
evaluated for the first time. The Galactic flux obtained are larger than those obtained 
previously in similar calculations. The non-annihilating scattering contributions of the 
propagated particles are introduced. The preliminary results are shown and discussed. 

1. Introduction 

A good knowledge of the Galactic production of (secondary) antiprotons 
and light antimatter nuclei, is a mandatory step on the way to the use of these (pri- 
mary) particles in the search for exotic astrophysical sources, like dark matter or 
primordial black holes [1]. These calculations all rely on: a) The empirical knowl- 
edge of the antiproton (antinucleon) production cross sections since no theory can 
predict them with the required accuracy to date; and b) On phenomenological 
models for the production of composite particles using these cross sections. In 
this work, the production of Galactic antimatter has been revisited on the basis of 
a new analysis of the experimental inclusive p production cross sections available, 
and of a reliable approach of the particle propagation with parameters tightly 
constrained by Cosmic Ray data. 

The atmospheric production of secondaries is a source of background for 
all embarked experiments (as well as a source of particles of high interest for 
other scientific fields). The knowledge of the corresponding flux in the Earth 
environment, is also a requirement for the interpretation of the satellite data 
which will be measured by the new generation of experiments. 
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2. Inclusive p production cross sections 

A selection of 654 data points of differential cross sections available from 
various pA —> pX experiments over a broad range of incident energies (typically 
10-450 GeV), has been analyzed consistently and fit using an analytical formula 
accounting for the cross section dependence on the reduced energy x* = E*/E^ nax , 
the transverse momentum (pt), and the nuclear mass of the target [2]. 

3. Coalescence models and antimatter nuclei 

The light antinuclei production cross sections were obtained by means of 
the microscopic coalescence model described in [3] . The simple coalescence model 
was also used for comparison. 

Recent nuclear collision data have accumulated experimental informations 
on the light nuclei and antinuclei production cross sections for various systems 
over a broad energy range. Some features particularly relevant to the present 
study are the quasi incident energy independence of the coalescence parameter 
for light nuclei production in p + A collisions, and the observed equality of matter 
and antimatter coalescence parameters [4]. A sample of available d data have 
been analyzed in this context, using the standard coalescence model. The results 
support the conclusions of the quoted reference. They will be reported later [8]. 

4. Particle production and propagation 
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Fig. 1. Simulation results for the atmospheric p (solid line) and d (dashed line) dux 
for the bins of latitude given on the figure. The d flux are multiplied by 2 I0 4 for 
convenience. 



3 



The secondary antinuclei production in the Earth atmosphere was gener- 
ated and propagated using a computer program developed recently and widely 
proven on various data samples (see [5] and refs therein). The resulting flux could 
be calculated at various coordinate positions covering the usual range of balloon 
and satellite experiments. The preliminary results for the p and d flux are shown 
on Fig. 1 for several bins of latitude between the equator and polar regions, at the 
AMS altitude (~ 380 km). The left and right panels correspond to the incoming 
and outgoing flux respectively. 

The p momentum spectra are peaked around 2 GeV per nucleon, which 
roughly corresponds to the velocity of the nucleon-nucleon center of mass system 
in which particles are produced with a maximum cross section at zero momentum, 
the incident energy times p production cross section distribution being peaked 
around 15 GeV. This spectrum is somewhat transformed by the cascade and 
transport in the atmosphere and in the Earth magnetic field, however. 



4-2. In the Galaxy 
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Fig. 2. Left: Calculated Galactic flux for d (solid line), t (dashed line), 5 He (dotted 
line), and 4 He (dash-dotted line) using production cross sections from the standard 
coalescence model and the LBM for the propagation. Similar results are obtained 
with the microscopic coalescence model. Right: Same as left for d (solid line) and 
with contribution of non-annihilating rescattering of the particles (dashed line), 
showing the dramatic increase of the flux at low energy expected from this effect. 
The dot-dashed and dashed lines correspond to the minimum and maximum flux 
respectively, constrained by the nuclear CR data [6], calculated in the Diffusion 
Model. 



The Galactic propagation of particles was performed using both the stan- 
dard Leaky Box Model (LBM), and the code developed by the LAPTH Annecy 
group, based on the Diffusion Model (DM) including convection, reacceleration 
and energy losses, with the propagation parameters constrained by the B/C ratio 
[6]. Fig. 2 shows the preliminary results. The left panel shows the calculated 
Galactic flux for a few light antimatter nuclei, d, t and 3 He, and 4 He. The flux 
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were calculated at TOA using the LBM, with solar modulation effects included 
($ = 500 MV). In the search of primordial antimatter, 4 He is the best candidate. 
It is thus of major importance to have reliable predictions for the Galactic sec- 
ondary flux of this particle, to which future measurements could be confronted. 
Fig. 2 right shows the same LBM d flux with the effects of non-annihilating scatter- 
ing of the propagated particles included. This latter contribution was calculated 
by using a non-annihilating inelastic cross section based on an analytical form 
for the pp — > pX differential cross section, renormalized using the total reaction 
cross section Od v j<y vv ratio. Note that this procedure is quite different than the 
usual ansatz. Details will be given in a forthcoming publication [8]. The DM 
calculations for d propagation are also shown on the figure (no rescattering effect 
included). They are in fair agreement with the LBM results, the latter being 
between the limits calculated by the former. The non-annihilating inelastic con- 
tribution in the DM approach, are being calculated. Note that the calculated d 
flux is found to be significantly larger, by a factor of about 2.5 than those obtained 
previously [7]. 

The ratio of the maxima of the calculated Galactic (~10~ 8 (m 2 .s.sr.GeV/n) -1 
and atmospheric d flux (see Figs. 1 and 2) is about 10. The atmospheric back- 
ground should thus not be a major source of experimental uncertainty since in 
addition, the knowledge of the particle momentum allows to discriminate atmo- 
spheric secondary particles from the Galactic ones. The d spectrum is found to 
peak around 4 GeV/n. This value is compatible with both a sal 6 GeV total en- 
ergy production threshold, and a peaking of the d producing incident protons in 
the 20 GeV range (see above). This is at variance however with the atmospheric 
d spectrum obtained, which is centered at a lower energy (2-3 GeV) . This feature 
is being investigated. 

Detailed calculations are in progress and will be reported later. 
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